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Abstract. We report the discovery of CFBDS J005910.90-01 1401.3 (hereafter CFBDS0059), the coolest brown dwarf iden- 
tified to date. We found CFBDS0059 using i' and z' images from the Canada-France-Hawaii Telescope (CFHT), and present 
optical and near-infrared photometry, Keck laser guide star adaptive optics imaging, and a complete near-infrared spectrum, 
from 1.0 to 2.2 yum. A side to side comparison of the near-infrared spectra of CFBDS0059 and ULAS J003402.77-005 206.7 
qq ' (hereafter ULAS0034), previously the coolest known brown dwarf, indicates that CFBDS0059 is ~ 50 ± 15K cooler. We 

, estimate a temperature of T e g ~ 620K and gravity of logg ~ 4.75. Evolutionary models translate these parameters into an age 

■ of 1-5 Gyr and a mass of 15 - 30 My,,,,. We estimate a photometric distance of ~13pc, which puts CFBDS0059 within easy 

reach of accurate parallax measurements. Its large proper motion suggests membership in the older population of the thin disk. 
The spectra of both CFBDS0059 and ULAS J0034 shows probable absorption by a wide ammonia band on the blue side of the 
//-band flux peak. If, as we expect, that feature deepens further for still lower effective temperatures, its appearance will become 
a natural breakpoint for the transition between the T spectral class and the new Y spectral type. CFBDS0059 and ULAS J0034 
would then be the first Y0 dwarfs. 
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(~ 100 000 K) to the coolest known brown dwarfs (pre- 
viou sly ULAS J003402 77-005206.7 (hereafter ULAS0034), 
>T8 Warren et al.l (12007b ). There remains however a sizeable 
temperature gap between the 600-700 K ULAS0034 and the 
~100 K giant planets of the Solar System. Many of the cur- 
rently known extrasolar planets populate this temperature in- 
terval, characterized by complex atmospheric physics: matter 
and radiation in these cold, dense, and turbulent atmospheres 
couple into a very dynamical mix, where molecules and dust 
form and dissipate. Current atmosphere models are rather un- 
certain in this unexplored temperature range and they will sig- 
nificantly benefit from observational constraints. Two major 
physical transitions are expected to occur between ~ 700 K 
and ~ 400 K and strongly alter the emergent near-infrared 
spectra ([Burrows et al.l l2003): NH3 becomes an abundant at- 
mospheric constituent and its near-infrared bands become ma- 
jor spectral features, and water clouds form and deplete H2O 
from the gas phase. The corresponding near-infrared spectral 
changes are likely to be sufficiently drast ic that the creatio n 
of a new spectra l type will be warranted ( Kirkpatrickl [2000). 
Kirkpatrick etaD d 19991) : iKirkpatrickl (|2000j) reserved the "Y" 



letter for the name of that putative new spectral type. 

To help fill this temperature gap, we conduct the Canad a 
France Brown Dwarf Su rvey (CFBDS dDelorme et al.ll2008l) ). 
which uses MegaCam ( Bo ulade et al.1 120031) /'and z' images 
to select very cool brown dwarfs (and high redshift quasars) 
on their very red i' - z' colour. We present here our coolest 
brown dwarf discovery to date, CFBDS J005910.90-01 1401.3 
(hereafter CFBDS0059), a >T8 dwarf with evidence for near- 
infrared NH3 absorption. Section 2 describes its discovery, and 
presents our follow-up observations: /', z', Y, J, H and Ks 
photometry and astrometry of CFBDS0059 and (as a refer- 
ence) ULAS0034, laser guide star adaptive optics imaging and 
a near infrared spectrum of the new brown dwarf. Section 3 
discusses the kinematics and the dynamical population mem- 
bership of CFBDS0059. Section 4 compares the spectrum of 
CFBDS0059 with those of Gl 570D (T7.5), 2MASS J0415-09 
(T8) and ULAS0034, and in the light of synthetic spectra uses 
that comparison to determine its effective temperature, gravity 
and metallicity. We also examine in that section the new spec- 
tral features which appear below 700K, in particular an NH3 
band, and discuss new spectral indices for spectral classifica- 
tion beyond T8. Finally, section 5 summarizes our findings and 
sketches our near-future plans. 

2. CFBDS 

Field ultracold brown dwarfs are intrinsically very faint, and 
as a result they can only be identified in sensitive wide-field 
imaging surveys. They are most easily detected in the near in- 
frared, and one could thus naively expect them to be most easily 
identified in that wavelength range. Brown dwarf spectra how- 
ever very much differ from a black body, and their considerable 
structure from deep absorption lines and bands produces broad- 
band pure near-infrared colours that loop back to the blue. At 
modest signal to noise ratios, those colours are not very dis- 
tinctive. Brown dwarfs are therefore more easily recognized by 
including at least one photometric band under 1 fan. At those 



shorter wavelengths their spectra have extremely steep spectral 
slopes, and the resulting very red i' -z' and z' -J colours easily 
stand out. 

As discussed in detail in lDelorme et all d2008l) . the CFBDS 
survey brown dwarf identification is a two-step process: 

- we first select brown dwarfs candidates on their red i' - z' 
colour in MegaCam images which cover several hundred 
square degrees. 

- /-band pointed observations of these candidates then dis- 
criminate actual brown dwarfs from artefacts, and astro- 
physical contaminants. 

The i' - z' selection takes advantage of the wide field of 
the MegaCam camera on the 3.6m CFHT telescope, and of 
the trove of deep observational material obtained with that in- 
strument. We rely on existing i' images from the Very Wide 
component of the Canada France Hawaii Telescope Legacy 
Survey (CFHTLS-VW) and, for different fields, on existing z' 
images from the Red sequence Cluster Survey 2 (RCS2). We 
then match those with either new z' or new i' exposures to ob- 
tain i' and z' pairs. The two parent surveys also obtain g' and 
r' images, which for the RCS2 survey are contemporaneous 
with the z' exposure. We don't use those as primary selection 
tools, but the contemporaneous exposures from the RCS2 sur- 
vey provide a welcome check that an apparently red source 
was not, instead, a variable brighter at all wavelenghts at the 
z' epochs. Since all our fields have low galactic extinction, the 
only other astrophysical point sources with a similarly red i' -z' 
are quasars at redshifts of z= a 6, which represen t the other mo- 



tivation of our program (e.g. Wil lott et al.l l2007). We discrimi- 
nate between quasars and brown dwarfs with 7- band photome 



try o btained on several 2 to 4m-class telescopes ( Del orme et al 



I2008L for details). These targeted follow-up observations also 
reject a number of remaining unflagged artefacts and they pro- 
vide a refined spectral type estimate, thanks to the much higher 
signal-to-noise ratio which we typically achieve on z' - / than 
on i' — z' ■ 



3. Observations 

We first identified CFBDS0059 as a brown dwarf candi- 
date when comparing a 360 s RCS2 z' exposure from 2005 
December 27 with a 500 seconds CFBDS 2' exposure from 
2006 August 31. CFBDS0059 is undetected in the i' image 
to i' AB = 25.2 (5<t), in spite of a strong z' detection (z' AB = 
21.93+0.05). The RCS2 survey obtains contemporaneous g', 
r' and z' images, and we checked the g' and r' exposures for 
a counterpart. These images, which were obtained within 50 
minutes of the z' observation, show no object at the position of 
CFBDS0059. This essentially excludes that the z' detection was 
due to a variable or slowly moving object with neutral colours. 
The i' AB -z' AB > 3.2 (5<x) lower limit was thus secure, and made 
CFBDS0059 a very strong candidate for follow-up. 



http://www-dapnia.cea.fr/Sap/Phys/Sap/Activites/Projets/Megacam 
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Fig. 1. K s -band finder chart for CFBDS0059. The field of view 
is 2 by 2 arcminutes. 



3.1. J -band photometry and near-infrared 
spectroscopy 

Our initial /-band imaging consists of five 20-second dithered 
exposures with the SOFI near-infrared camera on the ESO NTT 
telescope at La Silla on 2006 November 12. As discussed be- 
low, the photometric system of that instrument is non-standard. 
We used a modifie d version of the jitter utility within the ESO 
Eclipse package (iDevillardl 1 1 9971) to subtract the background 
and coadd the five exposures. We extracted photometry from 
the resulting imag e usin g PSF fitting within Source Extractor 
( Bertin & Arnouts ( 1996b and Bertin & Delorme, in prepara- 
tion) and obtain ,A/ ei , fl =18.1 1+0.06. The resulting zab - Jab = 
3.0 colour confirmed CFBDS0059 as a strong very late-T dwarf 
candidate a nd we triggere d //-band spectroscopic observations 
with NIRI dHodapp et alj|2003b at Gemini-North. Those were 
obtained in queue mode on 2007 July 30 and immediately con- 
firmed the very cool nature of CFBDS0059. We then requested 
J and /f-band spectroscopy with the same instrument, which 
was obtained on 2007 September 1 . 

All spectra were obtained through a ~ 077-wide slit, which 



produces a resolving power of 



A,) 



500. The object is 



dithered along the slit. The //-band spectrum is the sum of 16 
300-second integrations, while the Y + J and K band spec- 
tra each are the sum of 9 300-second integrations. Consecutive 
image pairs are pair subtracted, flat fielded using a median com- 
bined spectral flat and corrected for both spectral and spatial 
distortions. Spectra are extracted using a positive and a negative 
extraction box matched to the trace profile. A first wavelength 
calibration was obtained with argon lamp arc spectra taken at 
the end of the sequence, and the wavelength scale was then 
fine tuned to match the atmospheric OH-lines. Individual spec- 



tra extracted from image pairs were then normalized and me- 
dian combined into final spectra. Per-pixel S/N of 25, 40 and 5 
where achieved on the J, H and K-band peaks respectively. For 
all 3 wavelength settings the AO star HIP 105 12 was observed 
immediately after the science observations to calibrate the in- 
strumental spectral response and the telluric transmission. 

3.2. Additional near-infrared and optical Photometry 

The J filter of the SOFI camera on the NTT has a quite 
non-standard bandpass, for which the large colour correc- 
tions that res ult from the strong structrure in T dwarf 



spectra (e.g. Ste phens & Le ggett 2004) have not been fully 
characterized. To compare the spectral energy distribution 



(SEDs) of CFBDS0059 and ULAS0034 dWarren et al.ll2007l) 
we therefore prefered to obtain ad ditional near-infr ared 
wide-band photometry with WIRCam dPuget et alj 120041) on 
CFHT, whic h uses standard Mauna Kea Observatory in- 
frared filters ( Simons & Tokunaga 120021: iTokunaga et al . 2002; 
Tokunaga & Vaccal 120051 MKO system). The observations 
(2007 August 1st and 5th) used dithering patterns of ~ 60 arc- 
sec amplitude for total (respectively individual) exposure times 
of 300 (60), 150 (30), 300 (15) and 720 (20) seconds for the Y, 
J, H and K s bands. The skies were photometric and the seeing 
varied between 0.8 and 1.0". 

Table [Tj summarizes the magnitudes of CFBDS0059 and 
ULAS0034 in all available bands. The W I RCam photometry of 
ULAS0034 agrees with the lWarren et al.1 (120071) measurements 
within better than l cr for H band and within 1.5 cr for J. The 
Warren et al. (2007) K-band measurement used a K filter, while 
our WIRCam measurement uses the narrower and bluer Ks fil- 
ter. The 0. 15 mag difference betwe en these two observa t ions is 
approximately consistent with the IStephens & Leggettl (|2004) 
prediction for the effect of these different filter bandpasses 
at late-T spectral types. Similarly, the better short-wavelength 
quantum efficiency of the WIRCam detector can qualitatively 
explain the 0.2 mag nitude {2(f) d i screpa ncy between our Y 



photometry and the Wa rren et al.1 (120071) WFCam measure 



ment. The near-IR colours of the two brown dwarfs are similar, 
except H - Ks which is ~0.5 magnitude bluer for CFBDS0059 
than for ULAS0034. We will interpret the implications of this 
low Ks flux when we examine the near-infrared spectrum. 

3.3. Astrometry 

CFBDS0059 and ULAS0034 are serendipitously just 6.3 de- 
grees apart on the sky, and at similar photometric distances 
from Earth since they have similar spectral types and appar- 
ent magnitudes , and we initially entertained the idea that 
they might, perhaps, be part of a common moving group. The 
proper motion of CFBDS0059 however, measured between our 
2005 Megacam z' and 2007 WIRCam K s images, and uncor- 
rected for its parallactic motion, is fi a = +0.94 + 0.06"/yr, 
Hi = +0.18 ± 0.06"/yr (Table E). ULAS0034 on the other hand 
moves by u„ = -0 .12 ± 0.05"/yr and fi s = 0.35 ± 0.005"/yr 



dWarren et al] 120071) . The two proper motions are thus suffi- 
ciently different that the two brown dwarfs are clearly unre- 
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Table 1. Wide band i' AB and z' AB CFHT photometry measured with Megacam, M KO-system Y, J, H, and K s measured with 



WIRCam and associated colors. For comparison we include the War ren et al.l (12007b . Y, J, H, and K photometry of ULAS0034. 



Target 








Spectral band 








i' AB 


z > 


Y 


J 


H 


K s K 


CFBDS0059 


> 25.17±0.03 


21.93±0.05 


18.82±0.02 


18.06±0.03 


18.27±0.05 


18.63±0.05 


ULAS0034 


> 25.14±0.04 


22. 11 ±0.06 


19.12±0.03 


18.21±0.03 


18.49±0.05 


18.33±0.05 


Warren et al. (2007) 






18.90+0.10 


18.15+0.03 


18.49+0.04 


18.48±0.05 






Y — J 


J-H 


H-Ks 


H- K 




CFBDS0059 


3.11+0.05 


0.76±0.04 


-0.21±0.06 


-0.36±0.07 






ULAS0034 


2.99±0.07 


0.91+0.04 


-0.28±0.06 


0.16±0.07 






Warren et al. (2001) 




0.75+0.10 


-0.34+0.05 




0.01±0.07 





Table 2. 

2007) 



Astrometry of CFBDS0059. (Epoch:August 5th, 



RA(J2000) 
hh:mm:ss 



DEC(J2000) 
dd:mm:ss 



Us 



"/yr 



00:59:10.903 -01:14:01.13 0.94 ± 0.06 0.18 ±0.06 



lated. We checked for main sequence common proper motion 
companions to CFBDS0059, which would have provided wel- 
come age and metallicity constraints (e.g. iScholz et al.l 12003) 
but did not find any match within a 10 arc minute radius. 



3.4. Keck Laser Guide Star Adaptive Optics Imaging 

To search for binarity, we imaged CFBDS0059 on 
16 January 2008 UT using the laser guide s tar adaptive optics 



(LGS AO) system (Ivan Dam et all 120061: IWizinowich et al.l 
20061) of the 10-meter Keck II Telescope on Mauna Kea, 
Hawaii. Conditions were photometric with better than av- 
erage seeing. We used the facility IR camera NIRC2 with 
its narrow field-of-view camera, which produces an image 
scale of 9.963 + 0.011 mas/pixel dPravdo et al.ll2006l) and a 
10.2" x 10.2" field of view. The LGS provided the wavefront 
reference source for AO correction, with the exception of 
tip-tilt motion. The LGS brightness was equivalent to a 
V « 9.8 mag star, as measured by the flux incident on the AO 
wavefront sensor. Tip-tilt aberrations and quasi-static changes 
in the image of the LGS as seen by the wavefront sensor were 
measured contemporaneously with a second, lower-bandwidth 
wavefront sensor monitorin g the R = 14.6 m ag field star 
USNO-B1.0 0887-0010532 Monet et all 120031) . located 32" 



away from CFBDS0059. The sodium laser beam was pointed 
at the center of the NIRC2 field-of-view for all observations. 

We obtained a series of dithered images, offsetting the tele- 
scope by a few arcseconds, with a total integration time of 
1080s. We used the CHAs filter, which has a central wavelength 
of 1.592 fim and a width of 0.126 yum. This filter is positioned 
near the //-band flux peak emitted by late-T dwarfs. The im- 
ages were reduced in a standard fashion. We constructed flat 
fields from the differences of images of the telescope dome 
interior with and without continuum lamp illumination. Then 
we created a master sky frame from the median average of the 
bias-subtracted, flat-fielded images and subtracted it from the 
individual images. Images were registered and stacked to form 
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Fig. 2. Adaptive Optic detections limits as a function of the 
separation from CFBDS0059. Dashed lines show the equiva- 
lent mass and temperature for a given contrast. The first value 
assumes an age of 1 Gyr and the second an age of 5 Gyr. 



a final mosaic, with a full-width at half-maximum of 0.09" and 
a Strehl ratio of 0.05. No companions were clearly detected in 
a 6" x 6" region centered on CFBDS0059. 

We determined upper limits from the direct imaging by first 
smoothing the final mosaic with an analytical representation of 
the PSF's radial profile, modeled as the sum of multiple gaus- 
sians. We then measured the standard deviation in concentric 
annuli centered on the science target, normalized by the peak 
flux of the targets, and adopted 10cr as the flux ratio limits for 
any companions. These limits were verified with implantation 
of fake companions into the image using translated and scaled 
versions of the science target. 

Figure|2]presents the final upper li mits on any companio ns. 
We employed the "COND" models of iBaraffe etail d2003l) to 
convert the limits into companion masses, for assumed ages of 
1 and 5 Gyr and a photometric distance estimate of 13 pc. We 
assumed any cooler companions would have similar {CHAs - 
H) colors to CFBDS0059. 
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Fig. 3. Thin disk probability membership contours in proper 
motion space from the Besancon stellar population model. The 
contours are generated for synthetic stars with distances be- 
tween 10 and 20 pc, belonging to the thin disk (small dots) 
and the thick disk (small stars, with the density of the latter 
increased by a factor of 10 for display purposes). Less than 
one halo star would appear on the plot. Based on their mea- 
sured proper motion, the likelihood that CFBDS0059 (large 
open square) and ULAS0034 (large open triangle) belongs to 
the thin disk are 95% and > 99%. 



4. Kinematics 

We estimate a spectrophotometric distance for CFBDS0059 by 
adopting Mj - 17.5 ± 0.5, based on an approximate T9/Y0 
spectral type (discussed below) and o n an extrapolation o f the 
My versus spectral type relation of Kna pp et all (2004) be- 
yond T8 (2MASSJ0415 - 0935, hereafter 2M0415). Th is ex- 
trapolation is consistent with the IChabrier et al. ( 2000h mod- 
els which predict AJ ~ 1 .0 between brown dwarfs of T e ff ~ 
750K and 625K (like 2M0415 and CFBDS0059). The result- 
ing d = 13 + 5 pc has significant systematic uncertainties, be- 
cause spectral typing beyond T8 is just being defined, and es- 
pecially because the linear 1 magnitude/subtype decline seen 
at earlier subtypes may not continue beyond T8. The adaptive 
optics observations exclude any companion of similar luminos- 
ity beyond 1.2 AU, but CFBDS0059 could still of course be a 
tighter binary. Its small distance fortunately puts CFBDS0059 
within easy reach of modern parallax measurements; 



We use the Besancon stellar population model (IRobin et al 



2003) to generate synthetic stars between 10 pc and 20 pc for 
the thin (dots) and thick (star symbol) disk populations at the 
galactic position of CFBDS0059. Fig[3]shows their proper mo- 
tions together with those of CFBDS0059 and ULAS0034. The 
contour lines show the probability that an object with a given 
proper motion belongs to the thin disk rather than the thick disk 



(halo membership probabilities are negligible). CFBDS0059, 
at its probable distance, is well within the 95% probability thin 
disk membership region, and ULAS0034 is within the 99% 
probability region. In spite of its somewhat high proper mo- 
tion for an object beyond 10 pc, CFBDS0059 therefore most 
likely belongs to the thin disk. The mean age of the simulated 
stars in the region of the proper motion diagram occupied by 
CFBDS0059 is 4 Gyr, suggesting that it is an older member of 
the thin disk. That age is consistent with the 1 to 5 Gyr range 
derived below from comparison to COND evolutionary models 



(Bar affe et al.ll2003l) . As any kinematic age for an individual 



star, this determination has large error bars, but it suggests that 
CFBDS0059 might be older than ULAS0034. 

5. Spectral comparison and atmospheric 
parameters 

Fig|4]and Fig|5]present our spectrum of CFBD S0059, together 
with those of ULAS00 34 dWarren et al.ll2007l >T8). 2M0415 



dBureasser et al.1 120031 T8) and Gl 570D dBurgasser et al 
BOOOi T7.5), which successively were the coolest known 
brown dwarfs. Thanks to their earlier discovery, 2M0415 and 
Gl 570 have the best char acterized atmospheric parameters 
( Saumon et al. 20061 [2007b , and th ey provide the mos t solid 



baseline for a differential study. The Warren et al. ( 2007b spec 



trum of ULAS0034, kindly communicated by S. Leggett, was 
obtained with GNIRS on Gemini South and its jj = 500 reso- 
lution matches t hat of our NIRI spectrum of C FBDS0059. We 
downloaded the IBurgasser etail d2003i l2002h OSIRIS spec- 
tra of the two other brown dwarfs from the Ultracool Dwarf 
Catalog^, and degraded their original spectral resolution of 
^-1200 to match that of the GNIRS and NIRI spectra. 
Stronger telluric absorption from the lower altitude telescopes 
explains the wider blanked regions in the corresponding spec- 
tra, but doesn't measurably affect any comparison: as illustrated 
by CFBDS0059, late-T dwarfs have essentially negligible flux 
wherever telluric H2O absorption matters. Because the OSIRIS 
spectra do not cover the Y band, we complem ent t hem by lower 



resolu tion spectra from Geba lle et al.l (|2001) and Kna pp et al 
d2004 for A < 1.18/im. 



5.1. Atmospheric parameters 

Atmospheric parameters of ultracool dwarfs are ideally deter- 
mined from a combi n ation o f near and mid-IR information (e.g. 
Saumo n et al. l2006l I2007I) . but low resolution near-infrared 



spectra alone provide a useful proxy w hen mid-IR photometry 
and spectra are not (yet) available (e . g. IBurgasser et al.l l2006a; 



iLeggett et al.ll2007 ). Burgasser et al. ( 2006al) used a grid of so- 
lar metallicity cool brown dwarfs to calibrate two spectral ra- 
tios, H2O-/ and K/H, which respectively measure the strength 
of H2O absorption at ~ 1.15 /mi an d the flux ratio betwe en the 
K and H peaks, to T e ff and logg. IWarren et al.l d2007l) how- 
ever found that H2O-/ essentially saturates below T e ff = 750K, 
and therefore chose not to use this spectral index for spectral 
types later than T8. They demonstrate on the other hand that the 



2 http://www.iac.es/galeria/ege/catalogo_espectral/index.html 
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Fig. 4. 0.9 /an - 2.3 /an spectra of CFBDS0059 and the three other coolest brown dwarfs. The spectra are normalized to unit flux 
densities at their 1.27/an.peak, and vertically offset for clarity. The mai n T-dwarf spectral fe ature s are labeled. The tem peratures 
of Gl 570D and 2M0415 are from the careful spectroscopic analyses of lSaumon et al. (2006) and lSaumon et al. (2007). Those of 
CFBDS0059 and ULAS0034 are from our Wj versus J/K index (Fig.©. 



combination of the Kj J index with the width of /-band peak, 
parametrised by their Wj index, becomes a good T e ff and log g 
diagnostic at T e ff -900 K, and remains useful significantly be- 
low 750K. We adopt their method. 

Table [3] lists our measurement of these two indices for 



CFBDS0059, and Fig [6] compares them with the IWarren et al 



d2007l) measurements for Gl 570D, HD3651B, 2M0415 and 
ULAS0034. To derive T e ff and \ogg from the indices we 
use mode l indices from sola r -abundance BT-set tl atmospheric 
models dWarren et alj 120071: lAllard et all 120031 Allard et al. 
2008 in prep). The model with NH3 at chemical equilibrium 
abundance clearly produces too much absorpti on in the blue 
side of the i7-band, confirming the finding of ISaumon et alJ 
(2006, 2007) that non-equilibrium processes keep the NH3 par- 
tial pressure well below its equilibrium value. We then use 
models that keep the abundances of NH3 and N2 at a fixed value 
in all parts of the atmosphere where the reaction timescale ex- 
ceeds the mixing timescale, which typically occurs at the 600- 
800K temperature level. These "quenched" models agree much 
better with the observed band shape. 



As a first order correction for the remaining imperfec- 
tions of the theoretical spectra, the model indices are shifted 
into agreement with the measurements of 2M0415 at the 
[T eff =750K; log g=5.00 and rM/Hl=0 1 determined for that 
brown dwarf by ISaumon et al.1 (120071) . The T e ff=800K and 
log g=5.35 resulting from this calibration for Gl 570D are con- 
sistent with the T eff = 800-820K and logg=5. 1-5.25 derived by 
Saumo n et al. (2006) from a complete spectral analysis. For 
HD3651B, T eff =820-890K and ~ log g=5. 1-5.3 resulting from 
this calibration are roughly c onsistent with th e T e ff =780-840K 



and logg=5. 1-5.5 derived by lLiu et"aL I (l2007l) 

CFBDS0059 and ULAS0034 have very similar Wj indices, 
but the new brown dwarf has a significatively smaller K/J in- 
dex. Visual comparison of the two spectra (Fig.|4|i confirms that 
CFBDS0059 does have a weaker A'-band peak than any of the 3 
comparison cool br own dwarfs. As widely disc ussed in the re- 
cent lite rature (e.g. Liu et al. (2007), Fig . 3 in lBurgasser et alj 
d2006al) . or Fig. 3 in lLeggett et alj (12007b ). for a fixed metallic - 
ity a weaker /T-band peak is evidence for either a lower temper- 
ature or a higher gravity. The Wj index lifts this degeneracy: it 
indicates, again assuming identical chemical compositions for 
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Fig. 5. Y, J and //-bands spectra of the four cool brown dwarfs. The green vertical lines in the 1.15 /im - 1 .35 yum panel mark the 
two potassium doublets. The grey band in the 1.45 fim - 1.65 /im panel indicates the approximate limits of the new absorption 
bands, which we discuss in the text. 



the two brown dwarfs, that CFBDS0059 is cooler by 
and has a -0.15+0.1 higher log g than ULAS0034. 



-50+15 K 



As also discussed bv lWarren et alj d2007l) . the above uncer- 
tainties only reflect the random errors in the spectral indices. 
They are appropriate when comparing two very similar ob- 
jects, like CFBDS0059 and ULAS0034, since systematic errors 
then cancel out. They must on the other hand be increased to 
compute absolute effective temperatures and gravity: one then 
needs to account for the uncertainties on the 2M0415 parame- 
ters which anchor the Fig[6]grid (T~ ±25K and \ogg ~ +0.2; 



Sau mon et al.l d2007l) ). and for the uncertainties in the atmo- 



spheric models which may distort the grid between its anchor 
point [T eff =750K; logg=5.00 and [M/H]=0] and the ~600Kre- 
gion of interest here. We conservatively adopt T e ff=620+50K 
andlogg=4.75+0.3. 

This 2-parameter analysis obviously cannot determine all 
three main atmospher ic parameters CIW,, logg and metallic - 



ity). As discussed by IWarren et al.l d2007l) . it actually deter 



mines the temperature with no ambiguity but leaves a com- 
bination of [M/H] and logg undetermined, and they demon- 
strated that in the Wj versus J/K plot metallicity is degenerate 
with surface gravity, with A(logg) = -2A[M///]. CFBDS0059 
is thus definitely cooler than ULAS0034, but from Wj versus 
J/K diagram it can have either higher surface gravity or lower 
metallicity. This degeneracy affects the full JHK-band spec- 
trum, where any metallicity vs gravity difference is at most very 
subtle. It is however lifted by the sha pe of the T-band pea k 



(Figs. 3 of iBurgasser et al 



Legg ett et alj d2007)) 



d2006al) or 

since lower metallicity shifts the T-band flux density peak of 
submetallic brown dwarfs significantly blueward. Fig \5\ shows 



no such shift, and the two objects therefore have similar metal- 
licities. 

Fig. [7] overlays the observed CFBDS0059 spectrum with 
the synthetic spectrum for the closest point of the solar metal- 
licity atmospheric model grid. Except on the red side of the 
//-band, model and observations agree well, boosting our con- 
fidence in the derived atmospheric parameters. The main re- 
maining predictive shortcoming of the models is their overesti- 
mated absorption on the red side of the //-band peak. The prin- 
cipal opacity source in this region is the methane band centred 
at 1 .67 fim, for which comprehensive theoretical predictions are 
available, but only for transitions from the vibrational ground 
state (as will be discussed in detail in Homeier et al., in prepara- 
tion). To make up for the missing absorption from higher bands, 
which constitutes a significant fraction of the opacity at brown 
dwarf temperatures, a constant empirical correction factor was 
used. This correction must in turn lead to some overestimate of 
the CH4 absorption as we reach the lower end of the T dwarf 
temperature range. Another possible source of errors are uncer- 
tainties in the models' temperature profile. The BT-Settl mod- 
els self-consistently describe gravitational settling, convective 
turbulence, advection, condensation, coalescence and coagula- 
tion of condensates to pred ict the formation and vertical ex- 
tent of cloud layers dAllard et al.ll2003t iHelling et al1l2007l) . In 
late T dwarfs these clouds are predicted to reside deep in the 
optically thick part of the atmosphere. Their opacity is thus 
not directly visible in the spectrum, but it may still impact the 
thermal structure, and thus the relative abundance especially of 
temperature-sensitive species like CH4. 
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Fig. 6. Wj versus J/ K indices. The grid represents indices mea- 
sured on solar metallicity BT-settl model spectra, shifted into 
agreement of the [T eff =750K; logg=5.00 and [M/H]=0] with 
the measured indices of 2M0415. The error bars represent the 
lcr uncertainties of the measured spectral indices. 



Another (less serious) disagreement between the models 
and the observed spectra occurs in the Y band. The mod- 
els overestimate the flux on the blue side of the T-peak, and 
they imperfectly reproduce the general shape of the peak. As 
discussed below, the opacities in that band are dominated by 
pressure-broadened wings of the 0.77 fan K I line on the blue 
side and CH4 on the re d s ide. 

For iBaraffe et al.l (120031) evolutionary models, the 

4.45-5.05 determined above 



raffe et al.l (120031) 
T eff ^570-670K and logg 
translate into an age of 1-5 Gyr and a mass of 15My„ p (for 
1 Gyr) to 30M Jup (for 5 Gyr). The kinematics of CFBDS0059 
suggests that it belongs to an older population, and therefore 
slightly favour a higher mass and older age. 

5.2. Individual spectral features 

Fig. [5] zooms on the Y, J and //-band peaks of the four 
cool brown dwarf spectra. The published OSIRIS spectra of 
Gl 570D and 2M0415 do not cover the Y band, which instead 



is plot ted fr om the lowe r reso lution spectra of Geb alle et al 



d200ll) and iKnapp et al.l ([2004). For easier comparison, the 
CFBDS0059 and ULAS0034 F-band spectra are smoothed to 
that resolution. 

Direct comparison of the four spectra can be used to shed 
light on incipient new features and atmospheric chemistry. 
Features which are seen in both CFBDS0059 and ULAS0034 
are likely to be real even when their significance is modest in 
each object, and those which are absent or weaker in the two 
hotter brown dwarfs, can reasonably be assigned to low tem- 
perature molecules. Conversely, features which disapear in the 
two cooler objects trace higher temperature species. 



As discussed above, the F-band spectra of CFBDS0059 and 
ULAS0034 do not differ much. Given the strong sensitivity of 
that band to [M/H] that implies that the two objects have similar 
chemical compositions. The shape of the F-band peaks of these 
two coolest brown dwarfs on the other hand differ from that of 
Gl 570D and 2M0415, with the CFBDS0059 and ULAS0034 
peaks extending further into the blue. The dominant absorbers 
in the blue wings of the T-band peak is the pressure-broadene d 
wing of the 0.77 fim K I line (e.g. Burgasse r et al. 2006a). 
which must weaken as K I de pletes from the gas p hase under 
T e ft=~700K. As anticipated bv lLeggett et alj (120071) . the slope 
of the blue side of the T-band peak therefore shows good po- 
tential as an effective temperature diagnostics beyond spectral 
type T8. 

The strength of the J-band K I doublet is a goo d gravity 
estimator in ultracool dwarfs (e.g. lKnapp et alj|2004 ). because 
an increa s ed pre ssure at a fixed temperature favors KC1 over K 
(Lodders d 19991) 1 and consequently weakens atomic potassium 
features. At T eff ~ 750-800K the the /-band K I doublets re- 
main weakly visibl e and useful as a gravity proxy (Fig. 7 of 
IKnapp et all d2004 ). At T eff < 700K on the other hand, the 
K I doublets have completely vanished at the resolution of the 
current spectra (Fig. [5), even at the probably lower gravity of 
ULAS0034. Potassium is thus mostly converted to KC1 (or per- 
haps other compounds) in the relevant photospheric layers. 

The strongest new feature is wide absorption on the blue 
side of the //-band, at A < 1.565/im. It is conspicuous in 
CFBDS0059 and well detected in ULAS0034, and with hind- 
sight is weakly visible in the 2M0415 spectrum (Pig. [5j. It is 
however clearly stronger at T e ff < 700K. To visually emphasize 
this broad feature, we bin the spectra to R=~100 and overlay 
the four H-band spectra (Fig. [8] left panel). Absorption sets in 
at ~ 1.585jt/m and becomes deeper for A < 1.565yum. These 
wavelengths overlap with strong H2O and NH3 bands. Either 
molecule could a priori be responsible for the absorption. 

5.3. Near infrared ammonia signatures 

Those molecules unfortunately have imperfect opacity data, 
and the NH3 laboratory line lists in particular are incom- 
plete below 1.7yum. Comput ed ammonia opac ities are there- 
fore strictly lower limits. Legg ett et al.l (12007b compare syn- 
theti c spectra comp uted with and without NH3 opacity, using 
the llrwin et al.l (1 1999b line list for A < 1 .9/im, and find that am- 
monia absorption in cold brown dwarfs strongly depletes the 
blue wing of the H band (their Fig. 10). Similarly, Fig. 9 of 
Saum on et al. I d2000h plots synthetic //-band spectra with and 
without NH3 opacity, and find differences in two wavelength 
ranges: the NH3-rich model is significantly more absorbed for 
A < 1 .565^m and it has weaker but significant absorption in 
the [1.5725 - 1.585yum] range. 

Fig [8] right panel plots two BT-Settl models for [T e g = 
600K; logg = 4.75], without any near-infrared NH3 opacity, 
and with NH3 opacity for that molecule at its chemical equi- 
librium abundance. As discussed above the BT-Settl models 
do not reproduce the observed //-peak shape very well, and 
a quantitative comparison is thus difficult. The comparison of 
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BT-settl 600K; log g=4.75, quenched NH3 
BT-settl 650K; log g=4.75, quenched NH3 




A Gu.m) 



Fig. 7. Overlay of the CFBDS0059 spectrum with the solar metallicity [T eff =600-650K; log g=4.75] BT-settl synthetic spectrum. 
The two spectra are scaled to agree at their 1.27/vm flux peak. The "quenched NH3" models are chemical equilibrium models 
which enforce a constant abundance of ammonia in the cooler regions of the atmosphere. 



the two models nonetheless confirms the Saumon et al. (2000) 
conclusion that ammonia produces strong absorption below 
~ 1 .57yum and weaker residual out to 1 .595/im. These model 
predictions qualitatively match the behaviour seen in Fig[8] left 
panel. 

To emphasize the changes in brown dwarfs spectra when 
their effective temperature decreases from ~800 to ~600 K, 
we plot in Fig the ratio of the spectra of CFBDS0059 and 
Gl 570D. The signal to noise ratio of the resulting /f-band 
spectrum is too low for detailed analysis, and we therefore fo- 
cus on the Y, J and H flux peaks. To avoid confusion from 
changes in the temperature-sensitive methane bands, we also 
mostly ignore the parts of the spectrum affected by CH4 ab- 
sorption bands, hatched in dark and light grey for respectively 
stronger and weaker bands. Fig. [TOl shows the equivalent plot 
for ULAS0034, which is very similar. 

The //-band spectrum ratio prominently shows the new ab- 
sorption band, which outside the CH4 band closely match es the 
300 K NH3 transmission spectrum of llrwin et al. (1999). Both 
spectra are strongly absorbed between 1.49 and 1.52 jim and 



rebound from 1.52 to 1.57//m. Water absorption, by contrast, 
is a poor match to the features of spectrum ratio. The strongest 
water absorption (as computed from the HITRAN molecular 
database for a 600 K temperature) occurs below 1 .49/mi, at sig- 
nificantly bluer wavelengths than the CFBDS0059 absorption 
feature. 



Some weaker but still significant bands of the llrwin et al 



d 19991) laboratory ammonia spectrum occur in the J band. 
Those again match features of the CFBDS0059/G1 570D flux 
ratio, but that agreement is much less conclusive: water and 
ammonia absorptions overlap on the red side /-band peak, and 
CH4 absorption affects the blue side of the peak. A 1.25-1.27 
//m feature is seen on both flux ratios and on the ammonia trans- 
mission, and could be due to ammonia since it is clear of any 
strong H2O absorption band. The slight wavelength shift be- 
tween the laboratory and observed features however leaves that 
association uncertain. Detailed synthetic spectra based on fully 
reliable opacities will be needed to decide whether NH3 ab- 
sorption matters in the J band at the effective temperature of 
CFBDS0059. The main pattern in the T-band is a blue slope, 
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Fig. 8. Left: //-band spectrum of the four cool brown dwarfs binned to R=~ 100. The spectra are normalized at A = 1 .59/im. The 
integration intervals for the numerator and denominator of our proposed NH3 - H index are marked. Right: BT-Settl synthetic 
spectra for [T e ff = 600K; logg = 4.75] with and without near-infrared NH3 opacity; the NH3 abundance is at its chemical 
equilibrium value. 



which reflects the weaker pressure-broadened K I wing in the 
cooler brown dwarf. The weak 1 .03yum NH3 band is not seen. 

Ammonia is easily detected in mid-infrared SPITZER 
spectra for all spectral typ es cooler than ~T2 dRoellig et al. 
2004t ICushing et al.ll2006l) . though significantly weaker than 



initially expected because mixing from lower atmospheric 
levels reduces its abundanc e in the high atmo sphere below 
the local equilibrium value ( ISaumon et alJl2006l) . Weak near- 
infra red absorptio n by a mmonia has been tentatively detected 
by ISaumon et all d2000l) in the T7p dwarfs Gl 229B, but 
CFBDS0059 and ULAS0034 provide the first incontrovertible 
evidence for a strong near-infrared NH3 band in brown dwarf 
spectra. 



This conclusion contrasts with iWarren et al 
possible 



(120071) find- 
ammonia in 



ing possible but inconclusive evidence for 
UL AS0034. The main d ifference between the two analyses is 
that lWarren et alj (|2007) focused on a higher resolution search, 
at a necessarily lower signal to noise ratio, for individual NH3 
lines between 1.5 and 1.58//m. We instead looked for the global 



signature of the absorption band, which only becomes obvious 
when looking at the full //-band spectrum. 

5.4. Spectral type 

Table [3] lists for CFBDS0059 the spectral indices used by 
the spectral classification scheme of Burgasser et alJ (Eo0 6b). 
whic h refines the p r eviou s schemes of iGeballe et alJ (12002b 
and IBurgasser et al. (2002). These indices would imply a T8 
classification, identical to that of 2M0415. As discussed above 
however, the near-infrared spectrum of CFBDS0059 demon- 
strates that it is over 100 K cooler than 2M0415 and shows 
clearly different spectral features. Based on the new indices 
we present later, CFBDS0059 should be assigned a la ter spec- 
tral type. The almost identical Bur gasser et al. (2006b) indices 
of the two brown dwarfs instead reflect those indices measur- 
ing H2O and CH4 absorption bands which saturate and lose 
their effective temperature sensitivity at the T8 sp ectral type 
of 2M0415. Beyond T8 the Burg asser etaD d2006bl) classifica- 
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CFBDS 0059 and NH 3 transmission 



ULAS 0034 and NH 3 transmission 




CFBDS 0059 and H 2 transmission 



ULAS 0034 and H 2 transmission 




Fig. 9. Flux ratio between CFBDS0059 and Gl 570D (black), 
together with the l aboratory room te mperature transmission 



spectrum of NH3 dlrwin et al.l |1999) (red, top panel) and 



the 600 K H2O transmission spectrum computed from the 
HITRAN molecular database (red, bottom panel). The grey 
bands mark the parts of the spectrum affected by strong (dark 
grey) or moderate (light grey) CH4 absorption. 

tion scheme therefore needs to be extended, with new spectral 
indicators that do not saturate until significantly later spectral 
types. 

Fully defining this extension is beyond the scope of the 
present paper, since two known objects beyond T8 are not 
enough to explore spectral variability, but one can nonetheless 
start exploring. Since the main new feature is NH3 absorption 
in the blue wings of the //-band peak, we define a new NH3-Z/ 
index as 



NH 3 -H = 



(i) 



Its numerator and denominator range are plotted in Fig [8] 
The numerator integrates the flux within the main NH3 band 
and its denominator measures the bulk of the //-band peak 



Fig. 10. Flux ratio between ULAS0034 and Gl 570D (black). 
The overlays repeat those of Fig [9] 



(we note that in cooler objects, the denominator could be af- 
fected by some NH3 absorption; its integration boundaries 
might thus need to be refined after such objects have been 
discovered). We compute this index for Gl 570D, HD3651B, 
2M0415, ULAS0034 and CFBDS0059 (Table 0J, and find 
that it strongly decreases from G1570D to ULAS0034 and 
CFBDS0059 (which have very similar NH 3 -//). 

Over the limited effective temperature range spanned by 
G1570D, HD3651B, 2M0415, ULAS0034 and CFBDS0059, 
and as far as one can infer from just 5 examples, the NH3-H 
and Wy indices correlate strongly (Fig. fTTb. The numerator of 
W j is centered at wav elengths where both ammonia (Fig. 10 
of iLeggett et al.l (12007b ) and CH4 have significant opacity, and 
future modeling work should be able to establish whether the 
two indices probe the same molecule or not. Since the near- 
infrared spectra of ULAS0034 and CFBDS0059 differ signif- 
icantly more from that of the T8 2M0415 than the latter dif- 
fers from the T7.5 G1570D (as quantitatively demonstrated by 
Fig. ITTl. it is natural to assign a full spectral subtype to the in- 
terval between the two coolest brown dwarfs and 2M0415. By 
that reasoning, and if ULAS0034 and CFBDS0059 are consid- 
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Table 3. Measured s pectral classification indices for CFBDS0059 and ULAS 0034. The first six indices form the base of the 
Burgasser et alj d2006bl) classification sch eme, and the table in cludes the corresponding spectral classification on that scale. The 
Wj index is a recent addition proposed bv lWarren et al.l d2007l) . 
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UCIlUllllIlaLOI 


value 




spectral 


ValllG 




spectral 


indice 


wavelength 


wavelength 






type 






type 








CFBDS0059 




ULAS0034 


H 2 o-y 


1.14-1.165 


1.26-1.285 


0.029: 


t0.005 


T8 


0.012 


±0.007 


>T8 


ch 4 -7 


1.315-1.34 


1.26-1.285 


0.165: 


t0.005 


T8 


0.014 


±0.009 


>T8 


h 2 o-h 


1.48-1.52 


1.56-1.60 


0.119: 


t0.008 


>T8 


0.133 


±0.010 


>T8 


CH 4 -H 


1.635-1.675 


1.56-1.60 


0.084: 


t0.002 


T8 


0.096 


±0.006 


T8 


CH 4 -£ 


2.215-2.255 


2.08-2.12 


0.128: 


t0.037 


>T7 


0.091 


±0.015 


>T7 


K/J 


2.06-2.10 


1.25-1.29 


0.101: 


t0.002 




0.128 


±0.003 




w. 


1.18-1.23 


1.26-1.285 


0.257: 


t0.004 




0.262 


±0.006 





Table 4. NHy i / indic es for the five brown dwarfs discussed in this paper. The spectral classification either (1) follows 



Burgasser et al. (2006b) or (2) is derived from the NH3-H and W j indices, as discussed in the text 





G1570D 


HD3651B 


2M0415 


ULAS0034 


CFBDS0059 


NH 3 -/f 


0.672±0.008 


0.66±0.005 


0.625+0.003 


0.516±0.008 
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Spec. Type 


T7.5 (1) 
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Fig. 11. NH3-H index versus Wj index. The error bars represent 
the lcr uncertainties of the measured spectral indices. 



ered as T dwarfs, their spectral type should be T9, or perhaps 
slightly later. 

The T spectral class however is quite unlikely to remain 
the last spectral type, since for sufficiently low effective tem- 
peratures atmospheric models predict major changes in vis- 
ible and near-infrared brown dwarf spectra: NH3 bands are 
predicted to appear in, and eventually to dominate, the near- 
infrared spectrum, the strong pressure-broadened optical lines 
of Na I and K I are predicted to disappear as those atomic 
species get incorporated into molecules and solids, and water 
clouds are predicted to form and to largely deplete water from 



the gas phase dBurrows etal.ll2003t iKirkpatrickl koOS). since 
spectral classification, for mostly practical reasons, is tradition- 
ally based on optical and near-infrared spectra, such a major 
transition will justify the introduction of a new spectral type 



for which the Y letter ha s long been reserved ( Kirk patrick et al 
19991: lKirkpatrickl l2000). If the A ~ 1 .55 /mi NH 3 band keeps 



deepening as the effective temperature decreases further, and 
eventually becomes a major spectral feature, its appearance at 
T e s ^650 K will become a natural transition between the T 
and Y spectral classes. ULAS0034 and CFBDS0059 would 
then be the first Y dwarfs, and the prototypes for Y0 brown 
dwarfs, rather than T9. That decision will to some extent re- 
main a matter of convention, but it must in any case wait until 
larger numbers of similarly cool brown dwarfs can document 
spectral trends in finer detail, and preferably over a wider ef- 
fective temperature range. 

6. Summary and conclusions 

We have reported the discovery of CFBDS005 9, a very cool 



brown dwarf, discovered in the CFBDS survey dDelorme et al 



2008). Its effective temperature is ~50+15 K cooler than that 



of ULAS0034, most likely making it the coolest brown dwarf 
known at the present time. High spatial resolution imaging es- 
tablishes than CFBDS0059 has no similarly bright compan- 
ion beyond 0.09", and no companion with a contrast under 
3.5 magnitude beyond 0.3" (respectively 1.2 and 3.9 AU at 
the 13 pc photometric distance). Its kinematics suggest, with 
significant error bars, a ~4 Gyr age at which CFBDS0059 
would be a ~30Mj up brown dwarf. The atmospheric param- 
eters of CFBDS0059 however are compatible with any age 
from 5 Gyr down to 1 Gyr, for which its mass would be 
~15M j up . A trigonometric parallax measurement together with 
mid-infrared photometry and spectroscopy with SPITZER will 
significantly refine i ts physical parameters, as demonstrated by 
ISaumon et alJ (J2007) for slightly warmer brown dwarfs. 
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We assign absorption in the blue wing of the H-band peaks 
of both ULAS 0034 and CFBDS 0059 to an NH 3 band. If 
that assignment is confirmed, and if, as we expect, the band 
deepens at still lower effective temperatures, its development 
would naturally define the scale of the proposed Y spectral 
class. ULAS 0034 and CFBDS 0059 would then become the 
prototypes of the Y0 sub-class. 

The CFBDS survey has to date identified two brown dwarfs 
later than T8, CFBDS0 059 and ULA S 0034 ( which we iden 



tified independently of IWarren et al.l (120071) . iDelorme et al 



(2008)) in the analysis of approximately 40% of its final 
1000 square degree coverage. We therefore expect to find an- 
other few similarly cool objects, and hopefully one signifi- 
cantly cooler one. 

CFBDS0059 and ULAS0034 provide a peek into the at- 
mospheric physics for conditions that start approching those in 
giant planets, and the future discoveries that can be expected 
from CFBDS, ULAS, and Pan-STARRS will further close the 
remaining gap. They also bring into a sharper light the remain- 
ing imperfections of the atmopsheric models, and emphasize in 
particular the importance of more complete opacity data. Our 
analysis relies on a room temperature absorption NH3 spec- 
trum, but higher excitation bands than can be excited at 300 K 
must matter in T e s = 600 K brown dwarfs. The eventual iden- 
tification of ammonia absorption in the J band will also need 
complete opacity information for H2O and CH4 and full spec- 
tral synthesis, since the bands of the three molecules overlap in 
that spectral range. 

The spectral indices that define the T dwarf spectral class 
saturate below 700 K, and new ones will be needed at lower 
effective temperatures. We introduce one here, NH3-//, which 
measures the likely NH3 absorption in th e H band. Together 
with the Wj index of IWarren et al.l d2007l) and th e slope of the 



blue side of the F-band peak dLeggett et al.ll2007l) . it will hope- 
fully define a good effective temperature sequence. Metallicity 
and gravity diagnostics are less immediately apparent, but will 
need to be identified as well. 
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